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Abstract: Queueing problems are most common features not only in our daily-life situations such as bank 
counters, post offices, ticket booking centres, public transportation systems, but also in more technical 
environments such as in manufacturing, computer networking and telecommunications. For any queueing 
system average queue length, average system length, average waiting time in the queue and average 
waiting time in the system are the main observable performance characteristics. Control chart is a graph 
used to study how a process changes over time and it is also used to control ongoing processes. In this 
paper control limits are established to study the behavior of E k /M /I queueing model using the 
performance cliaracteristics. Numerical results are given to highlight its applications. 
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I. Introduction 

Queueing model with Erlang arrival has been discussed by Gross and Harris (1998) and several others. 
Control chart is a quality control technique evolved initially to monitor production processes. Montgomery 
(2005) proposed a number of applications of Shewhart control charts in assuring quality in manufacturing 
industries. Shore (2000) developed control chart for random queue length of M IM Is queueing model by 
considering the first three moments and also Shore (2006) developed Shewhart -like general control charts for 
G/G/S queueing system using skewness. Khaparde and Dhabe (2010) constructed the control chart for random 
queue length of M/M/l queueing model using method of weighted variance. Poongodi and Muthulakshmi 
(2012) analyzed number of customers in system of M/E k /1 queueing model using control chart technique. The 
Erlang arrival queueing system is applicable to many real-life situations: 

(i) In military recruits a recruit first lines up to have blood test at one station, an eye examination at the next 
station, mental test by a psychiatrist at the third, and is examined by a doctor for medical problems at the 
fourth and so on. Recruit needs to pass all phases before entering for selection. 

(ii) In a typical polling system, voters have to pass through many phases like showing identity cards, 
signing, getting ink mark etc., before voting. 

(iii) In an airline counter, passengers are expected to check in which consists of many phases, before entering 
into the plane. 

and so on. 

This motivated the author to study the construction of Shewhart control charts for number of customers 
in the queue, number of customers in the system, waiting time in the queue and waiting time in the system of 
E k /M/1 queueing model. 

II. E K / Mil MODEL DESCRD7TION 

Consider a queueing system in which the inter-arrival times follow k-Erlang distribution with meanl/X in which 
an arrival has to pass through k phases, each with a mean time llkX prior to entering the service. Service times 
follow an exponential distribution with meanl/ u. 

Let p n be the probability that there are n customers in E k /M/1 queueing system. Then 

Pn=T Pj (P) CD 

where p/ P) is the probability of completion of j phases and is given by 
ia n m 

klp^ roJ _ k a ^ k) 

where r 0 e (0,1) is the root of the characteristic equation 
u z k+1 - (kX+u) z + kX = 0 
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since p 0 = — taking p = A/u, equation (2) becomes 

k 

Pj (P) =p(l-r 0 )r 0 ^ 
Therefore equation (1), reduces to 
p„ = p (l-r 0 k ) (r 0 k ) n_1 
which is a geometric distribution. 

Control chart analysis of E k /M/1 queueing model is carried out in the following sections. 

III. Number Of Customers In The Queue 

Let L q denote the number of customers in the queue. 

The expected number of customers in the queue, E(L q ) is given by 

E(L q )= f>-l) p n 

= ^V(3) 
(l-r 0 ) 

and E(L q 2 ) = f>-l) 2 p n 

= pr 0 k (l + r 0 k ) 
(l-r 0 k ) 2 

Then variance of L q is 

Var(Lq) = E(L q 2 ) - (E(L q )) 2 
= pr 0 k (l+r 0 k -pr 0 k ) 
d-r 0 k ) 2 

Upper control limit (UCL), central line (CL) and lower control limit (LCL) of Shewhart control chart, under the 
assumption that the number of customers in the queue follows normal distribution, are given by 
UCL =E(L q )+3 JVar(L q ) , (5) 

CL = E (L q ) 



}" 

LCL = E (L q ) - 3 ^/VattLq) J 



The parameters of the control chart are obtained by using (3) and (4) in (5) as 



p r n k + 3 J p r n k (1 + r n k — p r n k ) 
UCL= — 5 VF 0 0 w 0 J 

CL=^!_ 
(1-ro ) 



TrT _ prp' -3Vpr 0 k (l+r 0 k -pr 0 k ") 



IV. Number Of Customers In The System 

Let L s denote the number of customers in the system (both in queue and in service). 
The expected number of customers in the system is given by 

E(L S )= f>p n 

= £np(l-r 0 k )(r 0 k ) n - 1 
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(1-ro ) 
and E(L S 2 ) = 2]n 2 p n 

n=l 

= P(l-ro k )En 2 (r 0 k ) n ^ 

= p(l+r 0 k ) 
d-r 0 k ) 2 

Then variance of L s is 

Var(L s )=E(L s 2 )-(E(L s )) 2 

= P(l+r 0 k ~P) (?) 
(l-r 0 k ) 2 

The parameters of Shewhart control chart, under the assumption that the number of customers in the system 
follows normal distribution, are given by 

UCL = E (L s ) + 3 ^/VarO^) , (g) 

CL = E (L s ) 

LCL = E (L s ) - 3 ^/VarCLs) 
Using (6) and (7) in (8), the parameters of the control chart are obtained as 



} 



ucl= p + 3Vp(l + r 0 k -p) 
(l-ro*) 



d-r 0 fc ) 
p-3Vp(l+r 0 k -p) 



V. WAITING TIME IN THE QUEUE 

Distribution of waiting time of a customer in the queue for the model under study is given by 

W q (t)= 1 -ro k e M(1_r 0 k)t , t>0. 

Then the pdf of waiting time is 

w q (t) = ur 0 k (l-r 0 k ) e" M(1 " r 0 k)t 

Let w q denote the waiting time of customers in the queue. 
The mean of w q is 

E(w q ) = u r 0 k (l-r 0 k )Jt e^ (1 ^° k) 'dt 



M(l-r 0 ) 

E(w q 2 ) = n r 0 k (1 - r 0 k ) J t 2 e^ (1 - r ° k) 'dt 
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!a 2 (l-r 0 k ) 2 
Then variance of w q is 

Var(w q )=E(w q 2 )-(E(w q )) 2 

= r 0 k (2-//) 
H 2 (l-r 0 k ) 2 



(10) 



The parameters of Shewhart control chart, under the assumption that the waiting time of customers in the queue 
follows normal distribution, are given by 

UCL =E(w q ) + 3 ^w,) , (n) 

CL =E(w q ) 

LCL = E (w q ) - 3 ^/VartWq) 
The parameters of the control chart are obtained by using (9) and (10) in (1 1) as 



} 



r 0 k +3Jr 0 k (2-r 0 k ) 

UCL= — v 0 , — - 

M(l-r 0 ) 



H(l-r 0 K ) 



r 0 k - 3Jr 0 k (2-r 0 k ) 
LCL= -° ^-JL^_ °_Z 

M(l-r 0 k ) 

VI. Waiting Time In The System 

Let w s denote the waiting time of customers in the system. 
The expected waiting time of customers in the system, E(w s ) is given by 
E (w s ) = E(L S )/ X 



Ml-r„*) (12) 
Also the variance of w s is given by 
Var(w s ) = Var(L s )/ X 2 

= P( 2 1+r ° k k - P 2 } (13) 
l 2 (l-r 0 k ) 2 

The parameters of Shewhart control chart, under the assumption that the waiting time of customers in the system 
follows normal distribution, are given by 

UCL = E (w s ) + 3 VVartwJ {U) 

CL = E(w s ) 

LCL = E (w s ) - 3 ^Var(w s ) 
The control chart parameters are obtained by using (12) and (13) in (14) as 
ucl= p + p(l + r 0 k -p) 
^d-r 0 k ) 

cl= ^nr- 

^(l-r G k ) 
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3^ P (l4 
X(l-r ( 



VII. Numerical Analysis 

Numerical analysis is carried out to analyze the performance of queueing system with reference to the 
parameters X, |a and k. As LCL values are negative for the selected values of the parameters, they are considered 
as zero and therefore not shown in the table as a separate column. 

Table gives the control chart parameters fornumber of customers in the queue, number of customers in the 
system, waiting time of customers in the queue and waiting time of customers in the system for certain selected 
values of X, u and k. 

Table - Control chart parameters for L q , L s , W q and W s 



Numerical values in the table reveal the following features: 
For all the performance measures of E k /M/1 queueing system 

(i) increase in the arrival rate X, increases the parameters and UCL for fixed values of n and k. 

(ii) increase in the service rate n, decreases the parametersCL and UCL for fixed values of X and k. 

(iii) increase in the number of phases k, decreases the parametersCL and UCL for fixed values of X and \x. 

VIII. Conclusion 

In this paper control chart technique is applied to analyze the system characteristics of E k /M/1 queueing 
system. These characteristics provide the arrivals to decide whether to join the system or not. In addition it gives 
the information on the waiting time. 
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